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Abstract 

We follow the evolution of the Ionization Potential (IP) for the paradigmatic quasi-one- 
dimensional trans-acetylene family of conjngated molecules, from short to long oligomers and to 
the infinite polymer trans-poly-acetylene (TPA). Our results for short oligomers are very close to 
experimental available data. We find that the IP varies with oligomer length and converges to the 
given value for TPA with a smooth, coupled inverse-length-exponential behavior. Our prediction is 
based on an “internally-consistent” scheme to adjust the exchange mixing parameter a of the PBEh 
hybrid density functional, so as to obtain a description of the electronic structure consistent with 
the quasiparticle approximation for the IP. This is achieved by demanding that the corresponding 
quasiparticle correction, in the GIT@PBEh approximation, vanishes for the IP when evaluated 
at PBEh(Q;*‘^). We find that is also system-dependent and converges with increasing oligomer 
length, allowing to capture the dependence of IP and other electronic properties. 
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I. INTRODUCTION 


Ionization potentials (IP) and electron affinities (EA) are fundamental electronic proper¬ 
ties of composite or complex systems and, notably in the past years when organic materials 
(molecular or polymeric) are being sought for applications in optoelectronic devices,-*^ much 
attention has been paid to this subject. The input from theoretical calculations is extremely 
relevant not only to help in gauging experimental data, but also in order to identify new 
directions for the optimal composition of actors in the building of a device. There has been 
thus a search for theoretical methods that can give us accuracy together with feasibility of 
calculations, spanning a wide range of both inorganic and organic molecular systems.-^- In 
particular, the class of linear or quasi-linear molecular systems -oligomers or polymers- offer 
a special work space that allows one to concentrate on length-dependence (just one relevant 
dimension) of the properties of interest: indeed, the dependence and evolution of the IP, EA 
and electronic gaps for with structural characteristics or compositions is a topic of intense 
studies.-^— For short oligomers in hxed geometries, as for small molecules, these properties 
can be obtained with high accuracy from high-level quantum-chemistry calculations that 
go beyond the mean-held approximation, serving as benchmarks for other computational 
electronic structure approaches.—'^ As the oligomer length increases, however, the compu¬ 
tational cost of such calculations quickly becomes prohibitive.— 

For such polyatomic systems, in particular for large molecules or extended materials, 
density-functional theory (DFT) has become the method of choice for a theoretical de¬ 
scription, analysis, or prediction of ground state electronic properties, stable or metastable 
atomic structures, vibrations, and structure-property relationships.— We recall that, de¬ 
spite the fact that DFT is a ground-state theory, certain excitations that can be expressed 
as differences of ground-state total energies are accessible. The IP and the EA are dehned 
as: 


( 1 ) 
( 2 ) 
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IP = 

EA = E^ - E^+^ , 



where ^ and E^~^^ are the total energies of the N-, (N-1)- and the (N+1 j-particle 

systems in the ground state. If E^ and E^^^ are computed for the same molecular geometry 
we obtain “vertical excitations”. The difference 

Egap = IP — EA (3) 

is the electronic gap of the system, also called the self-consistent or ASCF gap. Experimen¬ 
tally it is determined by direct and inverse photoemission, and should not be taken as the 
optical gap. 

In exact DFT the values of IP and EA from EqsiT]and[2]are also givenl^'i^ by the highest 
occupied Kohn-Sham (KS) levels of the N and N+1 electron systems, respectively. For 
approximate DFT, the Slater-Janak transition states, i.e. the highest occupied KS levels of 
the N-1/2 and of the N+1/2 electron systems, should provide an accurate estimate of the 
IP and EA energies: 


IP ^ 


N-l/2 


EA 


N+l/2 

^N+1 


( 4 ) 

( 5 ) 


The difference between the highest occupied KS levels of the N- and (N-l/2)- and of the 
(N+l/2)- and (N+l)-e\eciio\i systems reflects the self-interaction or localization error of 
the highest occupied KS orbitals of the N and (N+l)-e\eciio\i systems.— For approximate 
DFT the energies noted in Eqs. 0] and [5] should be taken.— >2^ 

The HOMO and LUMO (highest occupied and lowest unoccupied molecular orbital) levels 
of KS theory for the ground-state of a given A-electron system, that is, the e/j and 
energies, are however frequently used for the definitions of (the negative of) IP and EA, and 
the difference 


jpKS _ 

^gap ^N+1 


.N 


-N 


( 6 ) 


is usually termed the Kohn-Sham HOMO-LUMO gap. 

As said above, the use of beyond mean-field methods for large systems is still a challenging 
issue, and for this reason it is not known how the IP or the electronic gap develop as a 
function of polymer length. While Berger et al. showed,— through a “dielectric needle” 
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model for the polymer, that the polarizability per unit mer is inversely proportional to the 
polymer length, no such analytic dependence is known for the ionization potential. For 
approximate DFT functionals the IP from Eq. [H IPascF) is usually more accurate than 
from the plain HOMO energy —e^ because it is less affected by the self-interaction error.— 
However, going from short oligomers to more extended systems, different DFT functionals 
give rise to different length dependence of the IPascf, ranging from concave and straight 
to convex as a function of inverse length,— so the problem is still under discussion. 

To address this problem we will design a DFT functional consistent with many-body 
perturbation theory in the GW approach.— >21 GW has become the prime method for the 
computation of quasiparticle energies in solids as measured by direct or inverse photoe- 
missioii22i^— and is increasingly applied to organic systems^ii^ including polymers.— >2^^— 
The standard procedure is to apply a single iteration of the GW approach (GqWq) as a 
many-body perturbation to the results of a DFT or Hartree-Fock (HF) calculation. The 
single-particle wavefunctions of DFT or HF, the respective orbital energies and the resulting 
dielectric screening form the input to the GqWq calculation and therefore determine the be¬ 
havior of the screened Coulomb interaction Wq. The screening strength of W decreases with 
increase of the HOMO-LUMO gap. Thus local or semi-local DFT functionals that produce 
a too small gap compared to the real electronic gap would overestimate screening, whereas 
HF that produces a too large gap would underestimate it. As shown recently by Bruneval 
and Marques,— input orbitals and energies derived from hybrid functionals with a high frac¬ 
tion of exact-exchange yield GqWq IPs that agree well with experiment for small organic 
molecules, whereas for larger molecules the fraction of exact-exchange has to be consider¬ 
ably lower.— In other words, there is a well-known starting-point dependence of the GqWq 
approach,—^— and the best DFT starting point is usually also system dependent. It would 
therefore be desirable to iterate the GW approach towards self-consistency to eliminate the 
starting point dependence. Different schemes have been developed, either achieving self- 
consistency directly,—*^ or by the so-called ” quasi-particle self-consistency” of Schilfgaarde 
et al.,— which determines the variationally best non-interacting Green’s function Gq. The 
present work will follow a simpler and numerically more efficient approach. 

We here apply the GoHo@DFT approach to quasi-linear systems of increasing length. 
Trans-polyacetylene —(G 2 i? 2 )n— is the simplest conjugated material that already exhibits 
the alternating set of sp^-bonded carbon atoms, common to all conducting polymers, which 
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leads to tt- delocalization of the frontier molecular orbitals dictating the behavior of the 
electronic gap.— We will thus use the trans-acetylene (TA) family, from small oligomers 
(OTAs) to the inhnite polymer (TPA), as a model system to investigate the dependence 
of basic properties such as the IP and the HOMO-LUMO gap with localization length. 
To do that, we follow an approacb^>^ proposed recently: Building on the fact that the 
Kohn-Sham energy cat of the HOMO gives us the IP in exact DPT, we vary the amount of 
exact exchange in the Perdew-Burke-Ernzerhof hybrid functional^*^ (PBEh). We then pick 
that admixture a of exact exchange for which the KS-HOMO eigenvalue agrees with the 
quasiparticle energy from a GqWq calculation based on the same PBEh(Q;) starting point, 
denoted GoITo@PBEh(Q;*'^). The HOMO of PBEh is now consistent with the quasiparticle 
removal energy of GqWq and for this reason we call our scheme internally-consistent ic- 
PBEh. Monitoring for oligomers with increasing length then allows us to assess the 
length dependence of the IP, and gather information on the electronic screening. 

The remainder of the paper is organized as follows: In section m we provide a short 
overview of the basic concepts of the GqWq approximation and some technical aspects of the 
implementation; in section IHII we present our results, starting from the ground-state DFT 
calculations to obtain the geometrical models for the oligomers; we next present and discuss 
the internally-consistent model, applied to study the ionization potential of TA oligomers, 
with special focus on the length dependence. Finally, we draw our conclusions in section HVl 


II. THEORETICAL BACKGROUND 


In many-body perturbation theory the single particle excitation energies are the solutions 
of the quasiparticle equation 


+ Vext{r) VH{r)]'tpna{^) + / dr'S<^(r, r'; e® )?/;„^(r') = , 


(7) 


where Vext corresponds to the external potential created by the nuclei, vh is the Hartree 
potential, u is a state index and a the associated spin. The non-local complex self-energy 
operator S contains all electron-electron interaction effects beyond the Hartree mean held. 
In practice the self-energy needs to be approximated and we here adopt Hedin’s GW ap- 
proximation^*^ at the one-shot level 

, . i " 


■^GW, 


r,r ,ej = 


2ti 


de'G^(r,r',e + e')IFo(r,r',6')e*^'’ 


( 8 ) 
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where r] is an infinitesimal positive nnmber. Wo is the screened Coulomb interaction 


Wo{r,T',e) = j dr''e ^(r, r", e)n(r" — r'), (9) 

where n(r — r') = l/|r — r'| is the bare Coulomb interaction, and £“^(r,r";e) the inverse 
dielectric function. The latter can be written in terms of the polarizability 

£(r, r', e) = (5(r — r') — Jdr"v{r — r")Po{r", r'; e) (10) 

with 


Po(r, r'; e) = ^ J de'P^'^^G^ir, r'; e + e')G)((r', r; e'). 


( 11 ) 


Finally, Gq is calculated from the eigen-energies and wavefunctions of a preceding DFT or 
HF calculation 

^ ^ e - + IV sgniep - e^^)y ^ ’ 

Making the additional approximation that the quasiparticle wave functions equal the 
Kohn-Sham states, we can simplify Eq. [7] and write for the real part of the quasiparticle 


energies 


eZ = 4“ + » - ■K.c \;Pn,) = e'S + A". 


(13) 


where v^c is the exchange-correlation potential of the underlying DFT (or HF) calculation 
and (V'no-I {dy^(T)~'^xc It/'no-) the GqWq or quasiparticle correction. Equations [8] to 

Us] illustrate that and therefore the quasiparticle energies depend on the DFT functional 
used in the preceding calculation. 

As stated in the Introduction, in exact DFT the HOMO level of a finite system gives 
the IP and therefore the self-energy correction A^^qj^q is zero (for any other level no such 
statement holds). In standard approximations to the exchange-correlation functional the 
IP is typically not given accurately, because of the self-interaction error. Minimizing the 
absolute value of A^qj^q through optimization of a therefore implies that the self-energy 
correction to the HOMO level should be as small as possible. Alternatively, we could stay 
entirely within DFT and enforce the linearity of the DFT total energy with respect to the 
occupation of the HOMO stated to obtain a. However, this is not the scope of this paper 
and we will defer a discussion of the deviation of the straight line behavior and the internally 
consistent GW scheme to a forthcoming paper. 
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We emphasize that a is not related to a shift of the chemical potential (e^) that was 
originally proposed by Hedin,— who observed that if introduced in GqWq calculations it 
would model some effects of fully self-consistent GW calculations. The shift is also 
implemented in the GW space-time code^ and has negligible effects on the quasiparticle 
energies of semiconductors and insulators. This observation by Rieger et ah— is in line with 
the hndings by Pollehn et al,— who observe differences between GqWq, shifted GqWq and 
self-consistent GW only in the satellite spectrum of their Hubbard clusters and not in the 
quasiparticle peaks. 

To summarize this section, in the internally consistent GW scheme we explore the space of 
possible Go starting points spanned by the PBEh hybrid functional, and use the a parameter 
to traverse this space. In practice, we start from the same hybrid functional^ model of 
Perdew, Ernzerhof and Burke 

E,, = + (1 - 0 < a < 1 (14) 

where E^^ denotes the exact-exchange energy and and are the PBE exchange 

and correlation energy,— respectively. There the suggestion for a is 0.25, focusing on at¬ 
omization energies of a set of molecules. Here we follow a different rationale, and thus we 
perform a series of PBEh calculations for different values of a for the same molecule, and 
use the Kohn-Sham eigenvalues and orbitals as input for subsequent GqHo calculations. We 
hnd that [A^qj^qI can be minimized by just a few single shot GqWq calculations. 

III. RESULTS AND DISCUSSION 

In what follows we present Erst our results for TPA and OTAs obtained through different 
standard DFT functionals and HE, and discuss the convergence of (mean-held) electronic 
properties with conjugation length. Next we analyse the GoHo@DFT and GoHo@HF re¬ 
sults for differently sized OTAs, and proceed to the discussion of the internally-consistent 
procedure and the effects on the electronic structure in general. 

A. Starting-point calculations 

Our calculations are done for oligomers ranging from n=l (ethylene) to n=80 double 
bonds OTA(80) and for the inhnite TPA chain, over a single set of geometrical structures for 


all adopted functionals, so that we can evaluate the effect of each functional on the electronic 
properties independently from the effect on the structure. All calculations are performed 
using the FHI-aims code,— which has the advantage of including all electrons, a feature of 
basic relevance in our case as we will use the core-level energies explicitly for level alignment 
of different oligomers. Additionally, FHI-aims offers the possibility to calculate inhnite 
periodic as well as hnite systems with the same underlying approximations^ (e.g., basis 
sets, integration grid). FHI-aims is written with numerical atomic-centered orbital basis 
sets, organized in so-called “tiers” of basis sets, providing excellent convergence of density- 
functional based total energies even for complex structures, and sufficient convergence of 
GqWq results.— For geometric structure determination, we use the high accuracy tier 2 
basis set (Table 1 of Ref. 1^, 39 basis functions for C and 15 for H), and the calculations were 
carried out using the DFT functional of Perdew, Burke and Ernzerhof— (PBE) augmented 
by the Tkatchenko and Scheffler van der Waals scheme (vdW^‘^).— Unless otherwise stated, 
we select the tier 3 basis set (55 basis functions for C and 31 for H) to evaluate densities 
of states and quasiparticle energies. In the case of ASCF calculations we work within the 
spin-restricted, i.e. non-spin-polarized approximation to evaluate the total energies for the 
ionized systems. 

We hrst optimize the atomic coordinates for a TPA chain employing periodic boundary 
conditions, with the lattice constant along the chain direction hxed at the crystalline bulk 
value of c=2.457 A, as measured by X-ray scattering experiments.— To simulate a single 
isolated inhnite-polymer chain, the polymer backbone is placed in the (x,z) plane (the con¬ 
verged ground-state geometry of TPA is planar), and the lattice parameters perpendicular to 
the chain direction are set to a large value {a=h=25 A) in order to minimize the interaction 
between the chains in neighbor cells. A /s-point mesh of 1x1x10 is used in the optimization 
procedure. For these specihc settings, the resulting carbon-carbon bond distances are 1.362 
A for the double bond (C=C) and 1.423 A for the single bond (C—C), i.e., we obtain the 
expected dimerisation of the polymer backbone. The C—H bond length is 1.095 A and the 
C—C—C angle is found to be 123.8°. These structural parameters are quite similar to those 
found in previous theoretical studies using different functionals both within an oligomer 
approach^ or solid state calculations^'^ and also compares favorably with experimental 
results.— 
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TPA 





OTA n=8 


FIG. 1: (Color online) Schematic representation of the infinite polymer unit cell (TPA, top) and 
a finite model oligomer (OTA8, bottom), built by repetition of the unit cell. The hydrogen atoms 
added to saturate the oligomer chain and the resulting C-H distance are highlighted in red. 

The DFT equilibrium structure of the isolated TPA chain is then used as input to build a 
series of linear oligomeric chains OTA(n). Since the focus of the present work is the length 
dependence of the electronic properties of oligomers, we keep the relative atomic coordinates 
hxed at the inhnite chain result described above, and perform a further optimization only 
for the C—H distance of the end-cap CH 2 groups. The details of the geometry are not our 
main focus, as long as the geometry is consistent. We thus keep the atomic coordinates 
of the PBE-l-vdW^'^ optimization for each oligomer and apply different electronic structure 
approaches to these geometries. We hrst compare the following DFT functionals with HF; 
local-density approximation as parametrized by Perdew and Zunger (LDA-PZ),— PBE, 
and PBEO^. Then we perform GqWq calculations on top of these DFT functionals and on 
top of HF. We also apply our internally consistent scheme. The calculations are performed 
for the hnite oligomers with up to 30 double-bonds (n=30), which allows us to examine the 
length dependence of the frontier energy levels. 

Discussing hrst the results obtained with the standard PBE functional, we show in Fig. [2] 
the KS energy-level spectra obtained for a selected series of oligomers, and the density of 
states (DOS) for the inhnite polymer chain. To align the levels of all systems, oligomers and 
TPA, we use the average of the core Cis levels of each chain, that are then aligned at the 
value for the long oligomer OTA(50). As the chain length increases, we see the expected 
behavior of HOMO-LUMO gap closure, that converges to a small energy gap for the isolated 
TPA, in agreement with literature results.—*^ The main features of the continuum density 
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n=2 n=8 n=15 n=30 DOS (u.a) 


FIG. 2: Discrete DFT energy level spectra (obtained with the PBE functional) of trans-acetylene 
oligomers (OTAs), compared to the density of states (DOS) of the ID infinite polymer (TPA), 
calculated explicitly for the periodic model; here we use a gaussian broadening of 0.05 eV for the 
DOS. Spectra aligned at the average of the core levels (Cis) with those of OTA(50). 


of states of the polymer (i.e, the width of vr HOMO and LUMO bands, and the position of 
localized vr-states) starts to be visible for chains with ~n=15 double bonds, in agreement 
with previous theoretical estimates.- 

B. IP of the trans-acetylene oligomer series 

The scaling of physical properties of hnite conjugated oligomers as a function of 
chain length has been extensively studied and modeled in experimental and theoretical 
works,aiming to predict properties of polymeric materials using different extrap¬ 
olation models. We thus move now to the comparison with experimental results, and sum¬ 
marize in Table [I] the results for shorter oligomers, for which experimental data are avail- 
able.->^“— We hrst list the values coming directly from the negative of the HOMO eigenvalue 
(columns on the left) using the aforementioned different mean-held methods (KS and HF). 
Next we list the GqWq results for the corresponding starting point, and dually the results 
from the internally-consistent PBEh (columns on the right). We include also specihc litera¬ 
ture results obtained with the often-employed hybrid functional B3LYP.— The experimental 
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TABLE I: Ionization potential of acetylene oligomers, calculated at different theoretical levels: 
—Em negative of the KS (or HF) single-particle HOMO energy (left columns), our results for 
quasi-particle energies obtained through GoWo@DFT (right-columns), scGIT@HF, and from the 
internally-consistent procedure, see text (rightmost column). Experimental data included for com¬ 
parison at the center. All energies in eV. The mean absolute error for each functional as compared 
to the experimental values (n=l-4) is included in the last row. 


n 

DFT-KS DFT-GKS 

LDA PBE PBEO B3LYP^ 

HF 

Exp 

LDA 

GoWo 

PBE PBEO 

HF 

scGW icPBEh 

HF 

1 

6.85 

6.66 

7.77 

7.26 

10.10 

10.51-^ 

10.20 

10.25 

10.36 

10.70 

10.02 

10.44 

2 

5.95 

5.75 

6.68 

6.23 

8.60 

9.07-2^ 

8.65 

8.62 

8.83 

9.25 

8.47 

8.97 

3 

5.51 

5.31 

6.13 

5.69 

7.82 

8.29^^ 

7.78 

7.74 

8.00 

8.48 

7.65 

8.18 

4 

5.26 

5.04 

5.80 

5.36 

7.33 

7.79^ 

7.20 

7.17 

7.47 

8.01 

7.15 

7.69 

5 

5.09 

4.87 

5.57 

5.14 

7.01 

7.00^* 

6.84 

6.80 

7.10 

7.69 


7.36 

6 

4.97 

4.75 

5.41 

4.97 

6.78 

— 

6.56 

6.50 

6.81 

7.46 


7.12 

8 

4.59 

4.81 

5.19 

4.75 

6.48 

— 

6.17 

6.09 

6.44 

7.16 


6.79 

MAE 3.03 

3.23 

2.32 

— 

0.45 

— 

0.46 

0.47 

0.25 

0.19 

0.59 

0.10 


*Experimental (gas-phase) data available only for polyenes with terminal tert-butyl groups; in 
this Table, our calculated values for n=5 are also for butyl-terminated molecules. 

values are listed in the central column. 

Considering first the comparison between experiment and the LDA and PBE eigenvalues 
(—ej^) we see that, as expected, the gas-phase IP of all oligomers is strongly underestimated. 
The agreement with experiment is only slightly improved by the hybrid functionals (PBEO 
and B3LYP), while HF values are already very close. 

In Fig. |3](a,b) we select PBE, PBEO, and HF and now consider the evolution of the IP 
obtained from the difference in total energies (IPa 5 Cf)- In Fig. |3](c,d) we directly compare 
the PBE HOMO eigenvalue with IPascf for the PBE functional. We hrst note that IPa 5 cf 
from HF and DFT differ by a few electronvolts for small and medium sized molecules, but 
also that this difference tends to increase with oligomer length. The IPascf calculated with 
PBE (or PBEO) decreases quite fast with chain length. Indeed, we can see from Fig. [3] that 
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the slope of DFT-IPascf versus 1/1 increases with chain length, thus the value of the IP 
does not stabilize at longer chain lengths. A different trend is seen for HP, that is, the 
IPascf calculated with HP exhibits a decrease of the slope with growing oligomer length, 
a feature that can be seen more clearly following the inverse-length dependence. The fact 
that the slope in Fig. |3] increases for larger lengths in PBE and PBEO must be attributed 
to the semilocal part and not to the non-local exchange part, because it does not happen 
for HP. The difference between the DPT and HP IPascf reaches more than 1 eV in the 
inhnite chain limit. We see also in Fig. [3]^c,d) that the actual value of the negative of the 
PBE HOMO eigenvalue approaches the IPascf at the inhnite chain length limit, but the 
slope of the two curves are quite different. We observe the same behavior for PBEO (not 
shown here). 

Let us now proceed to the quasi-particle picture: while KS-HOMO levels of organic 
molecules are usually too high for local or semi-local DFT functionals, many-body corrections 
introduced perturbatively via GqWq calculations bring their values down, improving the 
description of IPs.— Our GqWq results for the IPs in Table |T] illustrate that the differences 
in the quasiparticle energies are indeed signihcantly smaller than the differences in the 
HOMO energy for the original DFT or HF values. The corrected values are all in much 
better agreement with the measured values, with a mean absolute error smaller than 0.5 
eV. However, contrary to the mean-held results, we now see an increasing deviation of the 
GqWo ip from measured values with increasing oligomer length. 

At this point it is illuminating to also inspect the self-energy correction A^qj^q to the KS 
HOMO level, as shown in Fig.|4]for PBE-based calculations. We observe that the self-energy 
correction decreases with chain length. This length dependence of A^^qj^q can be rational¬ 
ized in terms of a length-dependent change in the screening strength of the oligomer. Given 
the specihc vr-character of the frontier orbitals, the electron density of the KS HOMO state 
delocalizes over the backbone of the oligomers, and when the molecular length is progres¬ 
sively increased from OD ethylene towards quasi-lD oligomers, a substantial enhancement of 
the electronic screening is expected. The effect on the KS LUMO is similar for these systems 
as we will see, and thus this is reflected in the value of the electronic gap. Niehaus et al^ 
report similar conclusions for the band gap of ID polyacenes. Also for intrinsically differ¬ 
ent systems, sp^-bonded silicon nanocrystals, Delerue and co-workers^*^ observed in tight 
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Number of double bonds (n) Number of double bonds (n) 

2 4 6 10 15 20 25 20 10 6 4 3 2 



Number of double bonds (n) Number of double bonds (n) 


1 10 20 30 40 50 60 70 80 oo 20 10 6 4 3 2 



FIG. 3: (Color online) Top panel: Evolution of the first ionization potential of OTA series as a 
function of (a) chain length and (b) inverse chain length, up to n=30 double bonds, calculated 
through the ASCF approach with the different methods PBE, PBEO and HF. Bottom panel: Com¬ 
parison of the negative HOMO energy and the ASCF approach obtained with the PBE functional, 
in the (c) chain length and (d) inverse chain length representation, up to n=80 double bonds. The 
lines are just guides for the eye. 

binding GW calculations that the self-energy corrections to the DFT gap exhibit a smooth 
decreasing behavior with increasing (in that case 3D) nanocrystal size. These findings are 
in accordance with our results. 
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C. Internally-consistent mixing parameter 


We now move to the choice of the mixing parameter a to be inserted in the PBEh 
functional, Eq. [HI Fig. |5]shows the results for GqWq quasiparticle energies compared to the 
original PBEh(a) KS-HOMO energies, for three chosen OTAs n=2, 8 and 15. In this case, 
calculations are performed at the tier 2 basis set level. As a consistency check we compute 
the value for some selected oligomers using a larger basis set, namely tier 3, which allows 
for tightly converged orbital energies. The value is very stable with respect to the number 
of basis functions. Concerning the convergence behavior of both DFT eigenvalues and QP 
results we observe that the energies of the highest occupied states shift down by ^0.1 eV 
when going from tier 2 to tier 3 basis sets. These results indicate that tier 2 basis sets 
provide a good trade-off between accuracy and computational cost for the systems we study 
here. 

Figure [5] illustrates that the GqWo HOMO energy depends less on the a-parameter than 
the PBEh KS HOMO energy. The intersection between the PBEh and GqWq curves dehnes 
the internally-consistent fraction of EX namely and occurs at around a ~ 0.8. The 


Number of double bonds (n) 

1 5 10 15 20 25 30 



FIG. 4: Length (?) dependence of the self-energy correction on the HOMO energy A^qj^q of the 
OTA chains, defined as the difference between the quasi-particle and the KS energy; results for the 
PBE functional. The solid line is a fit with Ajj’qj^q = (1.6?“^/^ -|- 0.286)eV. 
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FIG. 5: Evolution of the Kohn-Sham and quasiparticle HOMO energy of trans-acetylene oligomers 
with increase of the mixing parameter a of PBEh. The parameter that satisfies the internal 
consistency criterion, indicated in the central panel, corresponds to the crossing-point between 
the curves calculated with PBEh(a) and GoIFo@PBEh(a). 


conjugated oligomers is thus much higher than the fraction included in most of 
standard hybrid functionals such as B3LYP (0.2), HSE or PBEO (0.25). 

Our zc-PBEh IP values are also included in Table U (zc-PBEh) for n=l to 8. ITe see 
that we obtain an improved description of the highest occupied state, yielding IPs in good 
agreement with gas phase reference data although no direct constraint is imposed in our 
scheme to fit experiment. Interestingly, the deviation of the zc-PBEh IPs from experiment 
remains approximately constant when increasing the oligomer length. This is shown in the 
inset of Fig. [6l We have also performed self-consistent GW (scGlE) for reference, using 
the scGW implementation in FHI-aims^>^. The results are included in Table [B The scGW 
IPs are consistently lower than experiment, GolTo@PBEO and zc-PBEh. This observation is 
consistent with recent benchmarks of scGW for molecules^‘^‘^“^‘^ and can be attributed 
to the pronounced deviation from the straight-line error (DSLE, also known as many-body 
self-interaction error) of scGlE— In GqWq, the DSLE can be reduced (or even eliminated) 
by an optimal starting point, which is why we here prefer to work with the ic-PBEh scheme. 

We now point out that the starting-point dependence of GqWq increases with system size, 
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FIG. 6: Chain length dependence of the first ionization potential of trans-acetylene oligomers 
obtained through Go Wo on different levels of mean-field methods: PBE (circles), PBEO (solid 
squares) and internally-consistent PBEh (stars); Hartree-Fock (empty squares). The lines are just 
guides for the eye. Included are also the experimental results (solid triangles) for small oligomers. 
Inset: difference between the calculated and experimental IP values for the small oligomers, same 
symbols as for the IP plots. 


reaching f^l.5 eV for / nm. This tells that the GqWo starting point is more important 
for longer or infinite chains. We thus plot in Fig. Hour calculated as a function of system 
size. We note that it varies slightly with chain length, ranging from a ~ 0.85 for the ethylene 
molecule (n=l) down to cc ~ 0.76 for the longest chain with n=30 double bonds. We ht 
with an exponential, as indicated in Fig. H which allows us to estimate the consistent 
fraction of exact exchange in PBEh as 75% for the case of an isolated TPA chain (we tested 
also htting with polynomials of (1//) but the errors are much larger). In contrast, Korzdorfer 
et al^ recently studied long-range hybrid functionals with an additional range separation 
between the long-range Coulomb potential and a short-ranged effective density functional. 
They chose to optimize the range-separation parameter, not an overall exchange-mixing 
parameter a, Ending that the range-separation parameter, optimized to satisfy the DFT 
analog of Koopmans’ theorem, strongly depends on the chain length and does not exhibit 
a saturation behavior for long polyene chains (n=25) whatever the nature of the starting 
functional. Thus, tuning the exact exchange parameter a as proposed by the internally- 
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consistent scheme^ is apparently a more adequate choice to predict the ionization potential 
for conjugated systems. As a last remark, we hnd that the optimal adjustment of the 
PBEh(Q;*'^) IP with length is obtained as IP(0 = a + h/l + {ce~^^)/l. A simpler inverse- 
length regression fails to reproduce the behavior at longer lengths as already discussed for 
energy gaps.— 

0.86 

0.84 

0.82 
b 

I 0.8 

0.78 

0.76 

0.74 

012345678 
chain length (nm) 

FIG. 7: Size dependence of the internally consistent a-parameter of the PBEh functional obtained 
for trans-acetylene oligomers. We see it decreases exponentially with chain length, as indicated by 
the dashed curve: the line is a fit with -|- 0.755. 

Coming to more specihc properties of these conjugated polymers, as pointed out before 
the KS LUMO shows vr-conjugated symmetry and localization properties similar to the 
KS HOMO - as a consequence we might expect that the behavior of the electronic gap also 
follows the same trend with tuning of a. We show in Fig. |8]the energy spectra obtained with 
our consistent procedure for intermediate length oligomers. We see that the HOMO-LUMO 
gaps obtained through PBEh(W’^) are close to the guasi-particle GqWq values. 

We now extend our analysis to the eigenvalues for a large energy window, and not only 
the KS HOMO itself, and in Fig. |9]we include the results just for the PBF, i.e. PBEh(a=0), 
and the PBEh(a®'^) starting points. In Fig. [9]^a) we show the A® energy corrections, for the 
same group of oligomers in Fig. |5l and in Fig. IHt^b) we show the complete spectra coming from 
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FIG. 8: Comparison between the energy level spectra calculated at the PBEh(a*'^) and the 
GoVFo@PBEh(a*'^) levels (with the internally-consistent mixing parameter) for the n=8, n=15 
and n=30 acetylene oligomers, close to the HOMO-LUMO gap energy window. 

the different methods, for the longer chain n=30 double bonds. Focusing on the eigenvalues 
close to the frontier orbitals, we first observe that the corrections to the first unoccupied 
PBEh(a!*^) states are really negligible, even for the n=15 double bond oligomer. Conversely, 
the corrections to the PBE eigenvalues are not only large (~1 eV) but also not constant as 
a function of energy. For the occupied states, the corrections to PBE are again large and 
not constant, but importantly they differ appreciably for localized and delocalized states. 
We recall that the peak in the density of states at eV below the KS HOMO in Fig. [H] 
derives from localized states, thus from that energy down to the valence band minimum 
we start to have mixing of states with different characteristics. The self-energy correction 
starting from PBE, for states in this energy window, is very different from the upper window 
values and not predictable, while the correction to PBEh(a*‘^) follows a smooth trend and is 
always less than ~1 eV; this is seen more clearly in Fig. [9](b). This analysis illustrates that 
PBEh(a!*^) is a more suitable starting point for GqWq than PBE, and thus also indicates 
that the PBEh(a*'^) spectrum is closer to experiment than PBE. 
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FIG. 9; Comparison of PBEh(a*'^) (ic-PBE) with PBE: (a) Deviation of the Kohn-Sham (KS) 
eigenvalues with respect to the GqWq quasi-particle spectrum plotted over the KS energies, for 
oligomer lengths n=2, 6 and 15 double bonds, (b) Spectra obtained for the longer oligomer n=30 
double bonds by the different methodologies. In (a) filled symbols indicate occupied states, empty 
symbols unoccupied states. 

IV. SUMMARY AND CONCLUSIONS 


In summary, we have presented and analyzed an “internally-consistent” [ic) parametriza- 
tion of the PBEh(a) functional that allows us to reproduce the electronic quasiparticle ener¬ 
gies normally obtained from GqWq calculations for the prototypical trans-acetylene family of 
conjugated systems. We show that the vertical ionization potential obtained with our opti¬ 
mized PBEh(Q;®'^)functional, that is, obtained through a non-empirical constraint, is always 
in much better agreement with available experimental values than if a simpler semilocal or 
standard hybrid functional (B3LYP, PBEO) is used. Furthermore, our internally consistent 
scheme also yields good KS LUMO energies that are consistent with GqWq, although this 
consistency is not a requirement in the construction of the scheme. We thus also hnd good 
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agreement for electronic gaps. As a last point, we show that many-body corrections to KS 
MO energies close to the frontier orbitals are also smaller than those obtained for standard 
functionals, which allows for sound prediction of the valence photoemission spectra. The 
dependence of the optimal internally consistent exchange mixing parameter a on the chain 
length is discussed and found to converge with increasing chain length. This is a signifi¬ 
cant computational advantage, as the controlled behavior of the a parameter should allow 
one to perform a single ic parametrization step for a certain class of systems and then use 
the corresponding PBEh(Q!*'^) functional for similar predictive simulations of the electronic 
structure of other, unknown systems of this class. 

Acknowledgment s 

The financial assistance of the Brazilian agencies CNPq and INCT-INEO are gratefully 
acknowledged. We also thank the action of CAPES (Brazil) and DAAD (Germany) for 
hnancial support under PhD grants. P.R. acknowledges the Academy of Finland through 
its Centres of Excellence Program (No. 251748 ). 


^ D. Cahen and A. Kahn, Adv. Mater. 15, 271 (2003). 

^ A. P. Kulkarni, C. J. Tonzola, A. Babel, and S. A. Jenekhe, Chem. Mater. 16, 4556 (2004). 

^ C.-G. Zhan, J. A. Nichols, and D. A. Dixon, J. Phys. Chem. A 107, 4184 (2003). 

Y. Zhao, N. E. Schultz, and D. G. Truhlar, J. Chem. Theory Comput. 2, 364 (2006). 

® G. Zhang and C. B. Musgrave, J. Phys. Chem. A 111, 1554 (2007). 

® I. Dabo, A. Ferretti, C.-H. Park, N. Poilvert, Y. Li, M. Cococcioni, and N. Marzari, 
Phys. Chem. Chem. Phys. 15, 685 (2013). 

^ P. Rademacher, K. Kowski, H. Hopf, D. Klein, O. Klein, and C. Suhrada, J. Mol. Struct. 567- 
568, 11 (2001). 

® A. Ruini, M. J. Caldas, G. Bussi, and E. Molinari, Phys. Rev. Lett. 88, 206403 (2002). 

® G. R. Hutchison, Y.-J. Zhao, B. Delley, A. J. Freeman, M. A. Ratner, and T. J. Marks, 
Phys. Rev. B 68, 035204 (2003). 


21 



C. Chi and G. Wegner, Macromol. Rapid Comm. 26, 1532 (2005). 

S. S. Zade, N. Zamoshchik, and M. Bendikov, Acc. Chem. Res. 44, 14 (2011). 

W.-F. Li, M. Andrzejak, and H. A. Witek, Phys. Status Solidi B 249, 306 (2012). 

A. Luzanov, J. Struct. Chem. 44, 681 (2003). 

M. Musia and R. J. Bartlett, Chem. Phys. Lett. 384, 210 (2004). 

M. Rohlfing, M. Tiago, and S. G. Louie, Synth. Met. 116, 101 (2001). 

R. B. Capaz and M. J. Caldas, Phys. Rev. B 67, 205205 (2003). 

A. Ferretti, A. Ruini, E. Molinari, and M. J. Caldas, Phys. Rev. Lett. 90, 086401 (2003). 

M. Levy, J. P. Perdew, and V. Sahni, Phys. Rev. A 30, 2745 (1984). 

C. -O. Almbladh and U. von Barth, Phys. Rev. B 31, 3231 (1985). 

W. Yang, A. J. Cohen, and P. Mori-Sanchez, J. Chem. Phys. 136, 204111 (2012). 

J. C. Slater (Academic Press, 1972), vol. 6 of Adv. Quantum Chem., pp. 1 - 92. 

D. A. Liberman, Phys. Rev. B 62, 6851 (2000). 

J. A. Berger, P. L. de Boeij, and R. van Leeuwen, J. Chem. Phys. 123, 174910 (2005). 

J. P. Perdew, R. G. Parr, M. Levy, and J. L. Balduz, Phys. Rev. Lett. 49, 1691 (1982). 

U. Salzner and A. Aydin, J. Chem. Theory Comput. 7, 2568 (2011). 

L. Hedin, Phys. Rev. 139, A796 (1965). 

L. Hedin and S. Lundqvist, Solid State Phys. 23, 1 (1969). 

W. G. Aulbur, L. Jonsson, and J. W. Wilkins, Solid State Phys. 54, 1 (2000). 

G. Onida, L. Reining, and A. Rubio, Rev. Mod. Phys. 74, 601 (2002). 

P. Rinke, A. Qteish, J. Neugebauer, and M. Scheffler, Phys. Status Solidi B 245, 929 (2008). 

X. Blase, C. Attaccalite, and V. Olevano, Phys. Rev. B 83, 115103 (2011). 

C. Faber, 1. Duchemin, T. Deutsch, C. Attaccalite, V. Olevano, and X. Blase, J. Mat. Sci. 47, 
7472 (2012). 

E. G. Ethridge, J. L. Ery, and M. Zaider, Phys. Rev. B 53, 3662 (1996). 

M. Rohlfing and S. G. Louie, Phys. Rev. Lett. 82, 1959 (1999). 

J.-W. van der Horst, P. A. Bobbert, M. A. J. Michels, G. Brocks, and P. J. Kelly, 
Phys. Rev. Lett. 83, 4413 (1999). 

J.-W. van der Horst, P. A. Bobbert, P. H. L. de Jong, M. A. J. Michels, G. Brocks, and P. J. 
Kelly, Phys. Rev. B 61, 15817 (2000). 

M. L. Tiago, M. Rohlfing, and S. G. Louie, Phys. Rev. B 70, 193204 (2004). 


22 


A. Ferretti, G. Mallia, L. Martin-Samos, G. Bussi, A. Ruini, B. Montanari, and N. M. Harrison, 
Phys. Rev. B 85, 235105 (2012). 

Y.-W. Ghang and B.-Y. Jin, J. Chem. Phys. 136, 024110 (2012). 

F. Bruneval and M. A. L. Marques, J. Chem. Theory Comput. 9, 324 (2013). 

T. Korzdorfer and N. Marom, Phys. Rev. B 86, 041110 (2012). 

T. Korzdorfer, R. M. Parrish, N. Marom, J. S. Sears, C. D. Sherrill, and J.-L. Bredas, 
Phys. Rev. B 86, 205110 (2012). 

N. Marom, F. Caruso, X. Ren, O. T. Hofmann, T. Korzdorfer, J. R. Chelikowsky, A. Rubio, 
M. Scheffler, and P. Rinke, Phys. Rev. B 86, 245127 (2012). 

P. Rinke, A. Qteish, J. Neugebauer, C. Freysoldt, and M. Scheffler, New. J. Phys. 7, 126 (2005). 
F. Fuchs, J. Furthmiiller, F. Bechstedt, M. Shishkin, and G. Kresse, Phys. Rev. B 76, 115109 
(2007). 

F. Caruso, P. Rinke, X. Ren, M. Scheffler, and A. Rubio, Phys. Rev. B 86, 081102(R) (2012). 
F. Caruso, P. Rinke, X. Ren, A. Rubio, and M. Scheffler, Phys. Rev. B 88, 075105 (2013). 

M. van Schilfgaarde, T. Kotani, and S. Faleev, Phys. Rev. Lett. 96, 226402 (2006). 

A. Moliton and R. C. Hiorns, Polymer Int. 53, 1397 (2004). 

V. Atalla, M. Yoon, F. Caruso, P. Rinke, and M. Scheffler, Phys. Rev. B 88, 165122 (2013). 

N. A. Richter, S. Sicolo, S. V. Levchenko, J. Sauer, and M. Scheffler, Phys. Rev. Lett. Ill, 
045502 (2013). 

J. P. Perdew, M. Ernzerhof, and K. Burke, J. Chem. Phys. 105, 9982 (1996). 

C. Adamo and V. Barone, J. Chem. Phys. 110, 6158 (1999). 

F. Bechstedt, Many-Body Approach to Electronic Excitations, vol. 181 of Springer Series in 
Solid-State Seiences (Springer, 2015), ISBN 978-3-662-44593-8. 

M. M. Rieger, L. Steinbeck, 1. White, H. Rojas, and R. Godby, Comput. Phys. Commun. 117, 
211 (1999). 

T. J. Pollehn, A. Schindlmayr, and R. W. Godby, J. Phys. Condens. Matter. 10, 1273 (1998). 
J. P. Perdew, K. Burke, and M. Ernzerhof, Phys. Rev. Lett. 77, 3865 (1996). 

V. Blum, R. Gehrke, F. Hanke, P. Havu, V. Havu, X. Ren, K. Reuter, and M. Scheffler, 
Comput. Phys. Commun. 180, 2175 (2009). 

X. Ren, P. Rinke, V. Blum, J. Wieferink, A. Tkatchenko, A. Sanfilippo, K. Reuter, and M. Schef¬ 
fler, New. J. Phys. 14, 053020 (2012). 


23 


V. Havu, V. Blum, P. Havu, and M. Scheffler, J. Comp. Phys. 228 , 8367 (2009). 

A. Tkatchenko and M. Scheffler, Phys. Rev. Lett. 102 , 073005 (2009). 

C. R. Fincher, C. E. Chen, A. J. Heeger, A. G. MacDiarmid, and J. B. Hastings, 
Phys. Rev. Lett. 48 , 100 (1982). 

P. A. Limacher, K. V. Mikkelsen, and H. P. Liithi, J. Chem. Phys. 130 , 194114 (2009). 

V. Lacivita, M. Rerat, R. Orlando, M. Ferrero, and R. Dovesi, J. Chem. Phys. 136 , 114101 

( 2012 ). 

S. Hirata, H. Torii, and M. Tasumi, Phys. Rev. B 57 , 11994 (1998). 

C. S. Yannoni and T. C. Clarke, Phys. Rev. Lett. 51, 1191 (1983). 

J. P. Perdew and A. Zunger, Phys. Rev. B 23, 5048 (1981). 

D. M. Ceperley and B. J. Alder, Phys. Rev. Lett. 45 , 566 (1980). 

S. Rohra, E. Engel, and A. Gorling, Phys. Rev. B 74 , 045119 (2006). 

™ S. Yang, P. Olishevski, and M. Kertesz, Synth. Met. 141, 171 (2004). 

S. S. Zade and M. Bendikov, Org. Lett. 8, 5243 (2006). 

J. Gierschner, J. Cornil, and H.-J. Egelhaaf, Adv. Mater. 19, 173 (2007). 

M. J. G. Peach, E. 1. Tellgren, P. Sack, T. Helgaker, and D. J. Tozer, J. Phys. Chem. A 111 , 
11930 (2007). 

J. Torras, J. Casanovas, and C. Aleman, J. Phys. Chem. A 116, 7571 (2012). 

M. Beez, G. Bieri, H. Bock, and E. Heilbronner, Helvetica Ghim. Acta 56, 1028 (1973). 

K. Kimura, Handbook of Hel photoelectron spectra of fundamental organic molecules: ioniza¬ 
tion energies, ah initio assignments, and valence electronic structure for 200 molecules (Japan 
Scientific Societies Press, 1981). 

T. Bally, S. Nitsche, K. Roth, and E. Haselbach, J. Am. Chem. Soc. 106, 3927 (1984). 

U. Salzner, J. B. Lagowski, P. G. Pickup, and R. A. Poirier, J. Comput. Chem. 18 , 1943 (1997). 
™ C. Rostgaard, K. W. Jacobsen, and K. S. Thygesen, Phys. Rev. B 81 , 085103 (2010). 

T. A. Niehaus, M. Rohlfing, E. Della Sala, A. Di Carlo, and T. Erauenheim, Phys. Rev. A 71 , 
022508 (2005). 

C. Delerue, M. Lannoo, and G. Allan, Phys. Rev. Lett. 84, 2457 (2000). 

G. Delerue, G. Allan, and M. Lannoo, Phys. Rev. Lett. 90 , 076803 (2003). 

F. Caruso, V. Atalla, X. Ren, A. Rubio, M. Scheffler, and P. Rinke, Phys. Rev. B 90 , 085141 
(2014). 


24 


M. Hellgren, F. Caruso, D. R. Rohr, X. Ren, A. Rubio, M. Scheffler, and P. Rinke, Phys. Rev. B 
91, 165110 (2015). 

M. Dauth, F. Caruso, S. Kiimmel, and P. Rinke, to be published. 

T. Korzdorfer, J. S. Sears, C. Sutton, and J.-L. Bredas, J. Chem. Phys. 135, 204107 (2011). 
H. Meier, U. Stalmach, and H. Kolshorn, Acta Polym. 48, 379 (1997). 

J. Rissler, Chem. Phys. Lett. 395, 92 (2004). 


25 


26 


